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The  acoustic  far-field  pressure  and  intensity  due  to  a  turbulent 
wake  is  calculated  using  Lighthill’s  free-turbulence  theory.  The 
radiation  intensity  is  plotted  as  a  function  of  the  Mach  number 
and  of  the  ratio  of  wake  wavelength  to  radiating  wavelength. 
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INFRASONIC  FAR-HELD  RADIATION  FROM  A  TURBULENT  WAKE 


INTRODUCTION 

The  problem  considered  in  this  short  report  is  the  calculation  of  the  far-field  intensity 
(or  pressure)  due  to  radiation  from  volume  turbulence.  The  specific  question  to  be  an¬ 
swered  is:  What  is  this  far-field  intensity  at  extremely  low  frequencies  (infrascund)? 

In  the  second  section  the  far-field  pressure  is  calculated  using  Lighthill’s  relation  for 
the  pressure  due  to  volume  turbulence  (1,2).  To  do  this  requires  certain  assumptions  about 
the  turbulent  velocity  field,  namely,  that  it  be  isotropic  and  have  the  form  of  an  outgoing 
spherical  wave.  The  basic  result  of  the  second  section  is  that  the  infrasonic  farfield  inten¬ 
sity  is  proportional  to  A/4 ( A/A)6  where  M  is  the  Mach  number,  A  is  the  turbulent  wave¬ 
length,  and  X  is  the  radiated  wavelength.  For  infrasound  (large  X)  this  factor  is  very  small. 
A  plot  is  given  showing  the  intensity  as  a  function  of  M  and  A/X. 

Finally,  the  last  section  contains  a  short  summary  of  both  assumptions  and  omissions. 


THEORY 

The  s|>ectral  pressure  p{ r)  at  the  field  point  r  (in  three  dimensions),  due  to  volume  tur¬ 
bulence  arising  from  a  turbulent  velocity  field  t>m(ro)  at  the  point  ro  (m  =  1,  2,  3),  is  given 
by  Lighthill  (1,2)  as 


Pn  C  c«fclr-*ol 

p(r)  =  4ff  |  aom  <Wm  ( r0 )^I ( 'o )  lr  _ yoy  dr0  (1) 


where 


30m  ‘  »io,„ 

differentiates  on  the  “zero”  coordinate.  The  Einstein  summation  convention  is  assumed 
here.  The  density  of  the  medium  is  po  and  k  =  2tt/X.  Further  manipulations  of  this 
volume  integral  yield  (1,2) 

Po  f  e/felf-roi 

Pir)  =  I  rfr0cm(ro)c„(r0)  jrT^j~  •  (2) 

In  the  far  field  (large  r) 


pP.’lr-rol  pikr  _ „ 

c  -  _  ^  £_  e-.!ir0vo»0o 

ir  —  r0 1  r 


(3) 
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where  r  •  ro  =  rr o  cos  Oq,  and  hence,  for  large  r, 

P(t)  ~  — ^Qmn  \4) 

where 

Qmn  -  4- J'^r0vm(ro)rfl(r0)e‘,7‘rocoftfo  .  (5) 

If  the  turbulence  is  assumed  to  be  isotropic, 

vm(r)  =  vm(r)  *  (6) 

and  behaves  like  a  spherical  wave  (with  wavenumber  K  which  could  have  a  positive  imagi- 
nar  part),  then 


ciKr 

Vm  (r)  =  Vm  • 


(7) 


Qmn  then  can  be  evaluated  as  (K  =  2ff/A) 

•2  n 


+oo  77  *  277 

Qmn  =  4^  J  dr°  J  d<?0  J  d0ofo2  sin  Oq  vm(r0)vn  (r0)e~ikrOcot0° 
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where  integral  tables  (3)  have  been  used  for  the  value  of  the  ro  integration.  Introducing 
tne  adiabatic  compressibility  *s,  defined  by 

A*  =  (Poc2)"1 

where  c  is  the  velocity  of  sound  in  water,  Eq.  (8)  becomes 
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Further  approximation  of  this  quadrupole  term  by  its  diagonal  values  yields  (6m„  is  the 
Kronecker  delta) 
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where  M  =  v/c  is  the  Mach  number  and  v  =  vm.  In  addition,  in  the  far  field,  the  term  in 
brackets  in  Eq.  (4)  can  be  written  as 


t)r  <)r 
OXm  dxn 


Using  the  approximation 


dr  dr 


and  the  far-field  intensity  /  as 


dxm 

dxn 

1 mn  =  &nn  ~  3, 
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For  very  low  frequencies  (infrasound,  X  large),  Eq.  (11)  can  be  written 

8  (fer)2JfJ“ 


(12) 


Hence,  in  addition  to  the  usual  spherical  spreading  loss,  there  is  a  dependence  on  the  Mach 
number  A/  raised  to  the  fourth  power,  and  to  the  ratio  of  turbulent  wavelength  to  radiated 
wavelength  raised  to  the  sixth  |x>wer.  There  is  also  a  factor  of  X  in  hr.  Figure  1  is  a  plot 
of  the  ratio  /(r)// o  versu*  A/X  from  Eq.  (12).  The  reference  level  is  a  comparison  source 
/o  having  an  output  of  0.64xl0“12  watts/cm2  (a  pressure  level  of  1  dyne/cm2).  The  quan¬ 
tity  I/1q  is  plotted  in  decibels  at  a  range  r  of  1  naut  mi  and  for  an  infrasonic  wavelength 
X  =  500  ft  (or  a  frequency  of  10  Hz).  The  intensity  at  a  new  wavelength  Xi  can  be  found 
from  the  relation  /j  =  /(X/Xj).  For  example,  doubling  the  wavelength  (Xj  =  2X,  frequency 
=  5  Hz)  decreases  each  of  the  10-Hz  curves  by  12  dB.  Halving  the  wavelength  (20  Hz)  in¬ 
creases  the  curves  by  the  same  amount. 


P-EMARKS 


There  are  several  assumptions  and  omissions  m  this  calculation  which  should  be  re¬ 
stated.  They  are 

a.  LighthilFs  turbulence  theory  for  free  turbulence  is  assumed  valid,  as  is  the  far- 
field  assumption.  Effects  due  to  boundary  layer  turbulence  have  been  omitted. 


